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Abstract. Using a quantum noise approach, we discuss the physics of both normal 
metal and superconducting single electron transistors (SET) coupled to mechanical 
resonators. Particular attention is paid to the regime where transport occurs via 
incoherent Cooper-pair tunneling (cither via the Joscphson quasiparticle (JQP) or 
double Josephson quasiparticle (DJQP) process). We show that, surprisingly, the 
back-action of tunneling Cooper pairs (or superconducting quasiparticles) can be used 
to significantly cool the oscillator. We also discuss the physical origin of negative 
damping effects in this system, and how they can lead to a regime of strong electro- 
mechanical feedback, where despite a weak SET - oscillator coupling, the motion of 
the oscillator strongly effects the tunneling of the Cooper pairs. We show that in this 
regime, the oscillator is characterized by an energy-dependent effective temperature. 
Finally, we discuss the strong analogy between back-action effects of incoherent Cooper- 
pair tunneling and ponderomotive effects in an optical cavity with a moveable mirror; 
in our case, tunneling Cooper pairs play the role of the cavity photons. 
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1. Introduction 

There has been considerable recent interest in studying the properties of mechanical 
oscillators coupled to quantum mesoscopic conductors. Such "quantum" nano- 
electromechanical systems (NEMS) are interesting because of their ability (in some 
cases) to perform quantum-limited position detection [TJ El El Gj and their potential 
to be used in quantum control applications In addition, they represent a new 
and interesting problem in the area of quantum dissipative systems, as the tunneling 
electrons (or quasiparticles or Cooper-pairs) in the conductor act as a non-equilibrium, 
non-gaussian dissipative bath for the mechanical oscillator. Theoretically, attention 
has largely focused on two systems: a normal metal single-electron transistor (SET) 
coupled to an oscillator [3 El El EEH El EEU U2] , and a tunnel junction (or quantum 
point contact) coupled to an oscillator [TSJ EH HI E] ; both systems have also been studied 
experimentally P3E1EII- ^ n these devices it has been predicted that for weak couplings, 
the back-action of the conductor essentially mimics the effects of an equilibrium thermal 
bath, with an effective temperature which is roughly proportional to the drain-source 
voltage in the conductor. 

In this paper, we turn to the properties of a quantum NEMS system which is 
considerably more complicated than either the normal- metal SET or tunnel junction 
systems. We consider a superconducting single electron transistor (SSET) coupled to a 
nanomechanical oscillator, focusing on regimes where transport in the SSET is via the 
incoherent tunneling of Cooper pairs (both the Josephson quasiparticle process (JQP), 
and the double Josephson quasiparticle process (DJQP)) [THl EI] . Such processes are 
attractive from a measurement point of view, as they have a high gain at a relatively low 
source-drain voltage. We also discuss the regime where transport through the SSET is 
due to the sequential tunneling of quasiparticles. Using a quantum-noise approach, we 
show that in each of these regimes the SSET provides the oscillator with a highly non- 
trivial effective environment. In particular, the effective temperature of the SSET can 
be much lower than the drain-source voltage in the SSET, meaning that considerable 
cooling can be accomplished; this is in stark contrast to the normal-metal SET and 
tunnel junction systems, where it is difficult to achieve cooling. In addition, the resonant 
nature of Cooper-pair tunneling leads to the possibility of unstable regimes characterized 
by negative dissipation; we discuss the strong electro-mechanical feedback that arises 
in these regimes and characterize the resulting stationary state of the oscillator. This 
state is highly non-thermal and can be described using an energy-dependent effective 
temperature. Finally, we point out that by using incoherent Cooper-pair tunneling, one 
can come extremely close to reaching the quantum limit on the displacement sensitivity; 
in contrast, this is not possible with a normal-metal SET. Note that an SSET plus 
oscillator system was also studied recently using an alternate technique by Blencowe, 
Imbers and Armour 18 . 

The plan of this paper is as follows. In section 2, we will review the effective bath 
description of a weakly-coupled NEMS, paying special attention to how back-action 
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damping arises both quantum mechanically and classically. In section 3, we introduce 
the SET plus oscillator system, and calculate the back-action in the sequential tunneling 
regime; our results here are applicable to a number of different systems, including 
a normal SET, a Coulomb-blockaded quantum dot, and a superconducting SET. In 
section 4, we calculate the back-action properties of incoherent Cooper-pair tunneling, 
and give a physical interpretation of negative damping effects that arise here. We point 
out the strong analogy between the superconducting SET system and an optical cavity 
with a moveable mirror, where forces arising from cavity photons can lead to cooling 
and instabilities [T[JJ 1201 EE] ; m our system, tunneling Cooper-pairs play the role of the 
photons. Finally, in section 5, we present results for the strong feedback regime that 
occurs as a result of these negative damping instabilities. 

2. Effective bath description and the origin of back-action damping 

At first glance, the theoretical problem presented by a NEMS system seems formidable: 
we need to understand how the back-action of our quantum conductor effects the 
oscillator, given that the conductor is typically in a far-from-equilibrium state, and 
typically does not have gaussian noise properties. Luckily, for many systems of interest 
the coupling between the oscillator and the conductor is sufficiently weak that the 
current in the oscillator responds linearly to the motion of the oscillator. As is discussed 
extensively in Ref. [221, the back-action effects of the conductor in this case can be 
directly related to the quantum noise properties of the uncoupled detector. We assume 
that the oscillator is coupled to the conductor via a term in the Hamiltonian which is 
linear in x: 

H mt = -Ax-F (I) 

Here, A is a dimensionless coupling constant; it needs to be sufficiently weak such that 
the conductor only responds linearly to the motion of the oscillator. The operator F 
describes the quantity in the conductor which couples to the oscillator; for the SSET, we 
will see that F is simply proportional to the charge on the central island. For simplicity, 
we will choose the oscillator's origin so that (F) — 0. F plays the role of a fluctuating 
back-action force; its unsymmetrized quantum noise spectrum at zero coupling to the 
detector, 

/oo A 
dte^(F(t)F(0)), (2) 
-oo 

will play a central role. 

Taking our oscillator to have mass m and natural frequency Q, and assuming that 
it is also coupled to an equilibrium Ohmic bath, one finds from perturbation theory in 
A that its motion is described by a classical Langevin equation [22]: 

mx(t) = -mn 2 x(t) - m lo x + 5f (t) + f avg (t) + 5f(t) (3) 

This classical Langevin equation involves two fluctuating forces, Sf (t) and Sf(t); it 
may be used to calculate the fluctuations in x{t) in terms of the spectral densities of the 
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Sfo and 5f fluctuations, as well as the oscillator's response to an external force. The 
answers obtained for both these quantities from Eq. ([3J) are in exact correspondence to 
what is found from a perturbative quantum mechanical calculation, performed to lowest 
non-vanishing order in A |22j . 

Consider the RHS of this Langevin equation. The second and third terms on the 
RHS describe the effects of the equilibrium bath on the oscillator: 7 is the damping rate 
induced by this bath, and 5f (t) is the corresponding fluctuating force. The fluctuation- 
dissipation theorem demands that the spectral density of the Sfo fluctuations obey: 

Ssf (w) = m^ohuj coth (hu/(2k B T )) (4) 

where To is the bath temperature. 70 and To have a simple interpretation: if the 
oscillator were not coupled to the conductor, its quality factor would be Qo = (the 
"intrinsic" quality factor of the oscillator) and its temperature would simply be T . 

The remaining terms on the RHS of Eq. Q describe the back-action of the 
conductor. Sf(t) is the fluctuating part of the back-action force; its spectral density is 
determined directly by the symmetrized quantum noise in the operator F: 

5*/H = + (5) 

In contrast, f avg (t) is the average value of the back-action force; it arises because the 
conductor (and hence (F)) changes in response to the motion of the oscillator. We can 
write this force as the sum of a conservative force which is in-phase with the oscillator's 
motion, and a damping force which is out-of-phase with the oscillator's motion: 

favgit) fdampit) -\- fconsit) (6^0 

/oo 
dt'-y(t - t')x(t') (6b) 
-oo 

/oo 
dt'a(t - t')x(t') (6 c) 

-oo 

Using standard quantum linear response relations, one has: 

X F (t) = -^9(t)([F(t),F(0)])- (7a) 
. . A2 f-hnX F (u)\ A 2 ( S F (u) - S F (—u)\ 

a(u) = A 2 (-Re\ F (uj)) (7c) 

Note that this description of the oscillator only requires A to be small enough that linear 
response is valid; there is no restriction on the size of the oscillator frequency Q. Note 
also that the fluctuating back-action force 5f(t) will not in general be related to the 
damping kernel 7(0;) by a fluctuation-dissipation relation, as our conductor is generally 
not in an equilibrium state. Nonetheless, we may use the fluctuation dissipation relation 
to define an effective temperature T e ff(u) at each frequency uj via: 

coth { , f \ = S -^l (8) 
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In the most general case, the oscillator responds to a variety of frequencies, and thus we 
cannot characterize the back-action by a single temperature; in this case, the oscillator 
will not be in a thermal state. However, if the oscillator has a sufficiently high quality 
factor, it only responds to forces in a narrow frequency band centered on uo = Q. If this 
frequency band is small enough, T e ff(u) will be constant over its width. In this case, 
the back-action will appear to the oscillator as being due to an equilibrium bath having 
a temperature T = T e ff(Q) and a damping rate 7(0). The spring constant k = mfl 2 of 
the oscillator will also be modified by the conservative part of the average back-action 
force; this will be given by: 

Ak = -a(ft) = A 2 (Re X F {Q)) (9) 

Including the effects of the equilibrium bath coupled to the oscillator, Eq. (J3J) tells 
us that the oscillator will be in a thermal state, characterized by a total damping rate 
7(fi) + 70 and a temperature T osc given rigorously by [22] : 

_ 7qT + jjtyTeff 

1 OSC — f s~\ \ \ L[J ) 

7o + 7(«) 

We again stress that the effective bath description presented above can hold even 
if the oscillator frequency Q is not small compared to the frequency scales of the 
conductor; all that is needed is that the "mechanical bandwidth" of the oscillator be 
small. Nonetheless, we are often interested in the case where I/O is long compared to 
the timescales relevant to the conductor; in this case, we can approximate 7(0;) ~ 7(0), 
a(tj) ~ ct(0) and Ssf(uj) — 5,5/(0). This corresponds to a damping force Fd am p = 
— 7777(0; = 0)x, a modified spring constant Ak = — a(u = 0), and a fluctuating back- 
action force with a white noise spectrum. In this case, Eq. (JEJ) reduces to: 

S F {0) A 2 S F {0 ) 
2d UJ S F (0)/h ~ 2m 7 

Before applying these ideas to the SET plus oscillator system, we comment on 
the origin of back-action damping, as there seems to have been some confusion on this 
point in the NEMS literature. Eq. (JT3j) makes it clear that there are two equivalent 
ways of thinking about damping. The first equality tells us that damping is simply the 
out-of-phase response of the average value of the back-action force F to the position of 
the oscillator. The second equality gives us a quantum picture for damping: damping 
results from the absorption of energy from the oscillator by the "bath" producing the 
force F (recall that positive-frequency noise corresponds to absorption of energy, while 
negative frequency noise corresponds to emission of energy |23j). 

There is another simple but useful classical way to understand both the origin of 
back-action damping and the spring-constant renormalization jTHlEOlE]- We simply 
need to use the fact that a) the average back-action force f aV g(t) is proportional to x, and 
b) the back-action force does not respond instantaneously to changes in the oscillator's 
motion. In the simple case where the oscillator is much slower than the source of the 



k B T eff EE _ = . V ' (11) 
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back-action, we can take the time delay to be given exactly by T resp \. We thus have: 

favgify favg\%(t ''"resp)] 

x(t Tresp) 

[x(t) - T resp x(t)} (12) 



df, 



avg 



dx 

df avg 



dx 

This immediately yields: 

Ak = (13a) 
dx 

dfavg /ion 

m7 = ~dx~ Tresp ' ' 

Thus, we can view the time delay in the response of the back-action force as being the 
origin of damping. Note that this argument tells us that Ak and 7 should always have 
the opposite sign, and that their ratio gives us a measure of r resp , the response time of 
the conductor's back-action: 

'resp — 

As we will see, the concept of a back-action response time will be useful in understanding 
the back-action properties of various different mesoscopic systems. 



3. Back-action of electron and quasiparticle sequential tunneling 

A SET consists of a metallic island with a large Coulomb charging energy Ec = e 2 / (2Cs) 
(Cs is the total capacitance of the island) coupled via tunnel junctions to both a source 
and a drain metal electrode; in the case of a SSET, the island and both the electrodes 
are superconducting. The charging-energy term in the Hamiltonian is given by: 

H c = E c {h-M)\ (15) 

where n is the charge on the SET island, and M = C g V g /e is the dimensionless electron 
number associated with a gate voltage V g which is coupled to the island via a capacitance 
C g . In addition, a voltage Vsd is applied between source and drain which drives the 
tunneling of electrons across the SET. 

To make a NEMS device involving a SET, a mechanical oscillator (coated with 
metal) is placed in proximity to the island of the SET such that there is a capacitive 
coupling C osc between the charge on the SET island and the potential of the oscillator; 
C osc depends on the island-oscillator distance. The oscillator is then voltage biased, 
with the result that changes in its displacement modify the electrostatic potential of the 
SET island (i.e. the parameter M becomes x dependent). For small displacements x, 

I A more exact equation would be to say f a vg(t) — %^ I- 00 ^ ex P[ — — t')/T resp ]x(t r ). For 
QT resp <C 1, this reduces to Eq. I|12|) above. 
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the result is a linear-in-x coupling term in the Hamiltonian 



int 



— Anx 



(16) 



where A is the effective coupling strength: 




osc 



1 

d 



(17) 



Here, V osc is the voltage applied to the oscillator, and d is the oscillator-island spacing 
when the oscillator is in its equilibrium position; for more details, see, e.g., Ref. [TO] . 
Comparing against Eq. ^ we see that Ah plays the role of a fluctuating back-action 
force on the oscillator, and that it is the quantum noise spectrum of the charge h that 
will determine the nature of this back-action. This conclusion is true regardless of where 
the SET is operated, or whether it is superconducting or normal. 

Calculating the back-action thus reduces to a problem of calculating the quantum 
noise spectrum of charge fluctuations in the SET. For various operating points of a 
SET/SSET, this has been done |21 ESI EH1 EH EB1 E3] ■ The simplest operating regime 
of a SET is that of sequential tunneling. Current flow is due to electrons (or, in the case 
of a SSET, quasiparticles) incoherently tunneling on and off the SET island via real, 
energy conserving transitions; because of the large charging energy Ec, only two charge 
states of the island are involved in transport. The quantum charge noise in this regime 
was calculated in [23] by using a qubit as a theoretical quantum noise spectrometer. 
In this approach, one first models the dynamics of a two-level system (TLS) weakly 
coupled to the SET island charge, and calculates an equation of motion for a reduced 
density matrix tracking both the TLS state and the SET island charge; for sequential 
tunneling, this is done to lowest non-vanishing order in the tunneling in the SET. In 
the limit of weak SET-TLS coupling, one can calculate both the stationary state of the 
TLS, as well as its relaxation rate. These may then be used to directly extract the 
finite frequency quantum noise spectrum of n at the splitting frequency of the TLS. By 
varying this splitting frequency, one can extract the entire quantum noise spectrum of 
h. 

Using the above approach, one finds that the low-frequency SET charge noise in the 
sequential tunneling regime is simply telegraph noise, while the frequency-asymmetry 
in the noise is determined by how the tunneling rates depend on energy [2E1 EH] • To 
be definite, consider a sequential tunneling process involving the island charge states 
n = and n — 1, and consider the zero temperature limit. There are then two rates of 
interest: r + , the rate at which electrons (or quasiparticles) hop onto the island from the 

§ Note there is an additional term in the SET plus oscillator hamiltonian which is proportional to 
V D 2 sc x 2 , and which thus contributes a V^-dependent shift in the oscillator's frequency. This term has 
nothing to do with back-action: it is independent of the SET island charge n. As such, we do not 
discuss it in what follows. 
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left junction, and T_, the rate at which electrons hop off the island through the right 
junction. The average current is given by I = er + r_/T s (where T s = T + + r_), and 
the average island is (n) = eT + /TY,. It will also be important to know how these rates 
change if we either increase or reduce the amount of energy driving the given tunnel 
event; let E represent this additional added energy Using Eqs. (|76jl and (fTTjl and the 
results of Ref. [23J, one obtains the following simple expressions for the back-action in 
the limit of small oscillator frequency (i.e. Q <C T^): 

r + r_ 
d E (r + r_) 



k B T ef f = - ^ - ; (18a) 



-7 = A 2 ^±^ (186) 
5k = - (mi) x T s (18c) 

Here, each rate should be evaluated at E — 0. These expressions are equally valid for 
sequential tunneling of electrons in a SET, electrons in a Coulomb-blockaded quantum 
dot, or BCS quasiparticles in a superconducting SET ||: the only difference is in the 
form of the tunnel rates. Note that the bath response time r resp defined in Eq. (jUj) 
is simply given by r resp = 1/(Te), as could have been expected: this is the timescale 
which characterizes charge relaxation on the SET island. 

For a normal-metal SET, in the case where both junctions have a conductance 
ge 2 /h and equal capacitances, the tunnel rates are given by 

r±(£) = 9{E± + E) [1 + n B (E ± + E)) , (19) 

where rig is the Bose-Einstein distribution function, and E± is the gain of electrostatic 
energy associated with each tunneling process: 

E + = eV ds /2 - 2E c {\/2 - J\f) (20a) 
E_ = eV ds /2 + 2E c (l/2 - AT) (206) 

As a result, at zero physical temperature in the SET, the effective temperature and 
back-action damping are given by: 

k B T eff = poPieV ds (21a) 
m 1 =A 2 -—^— (21b) 

g(eV ds y 

where po = T^/Ty, is the probability of having n = 0, p\ = 1 — po is the probability of 
having n — 1. For a fixed V^, the damping is constant, while the effective temperature 
reaches its maximum at the point of maximum current, M = 1/2. Note these 

results are in agreement with Refs. [01 H2] , which treats the SET-oscillator system using 
a generalized master equation. 

|| Note that there is a difference between a normal-metal SET and a Coulomb blockaded quantum 
dot: the former system has a vanishing single-particle level spacing, while the latter system has a large 
single-particle level spacing. As a result, the tunneling rates and back-action are quite different. 
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Turning to the case of a superconducting SET (i.e. SSET) with identical junctions, 
the quasiparticle tunneling rates are now given by: 

r± = Iss [(E± + E)/e] [1 + n B {E ± + E)\ (22) 

Here, IssiY) is the usual I-V characteristic of a SIS junction [2H]: there is no current 
until V = 2A/e, after which there is an abrupt rise. In the usual case where there is 
some inelastic scattering in the superconductor, this discontinuity in the I-V curve is 
smeared out over a voltage hj (e ■ tq P ) <C A/e, where tq P is the quasiparticle lifetime in 
the superconductor [SOI- If Note that the threshold voltage for quasiparticle sequential 
tunneling is Vds = 4A, as one needs to create four quasiparticles to transfer a charge 
e from left to right. To compare against the normal SET case, we further specialize to 
the case where M = 1/2; this gives a maximum current. We now have for both normal 
and superconducting SET's that T + = T_ = T, and: 

k B T eff = (23a) 



For drain source voltages much larger than 4 A/e, the tunnel rates for the quasiparticles 
will be almost identical to that for normal electrons, and the back-action will be almost 
identical to that in a normal SET at a similar voltage. However, for eVns ~ 4 A, the 
sharp rise of the SIS current-voltage characteristic will lead to stronger back-action in 
the superconducting case. One finds that c^T will be larger than in the normal state 
case by a large factor Atqp/Ii. One thus has: 

k B T eff ~ — (24a) 

TQP 

mi ~A 2 ^ (246) 

Thus, near threshold, quasiparticle sequential tunneling can lead to very low effective 
temperatures and high damping rates. This behaviour is entirely due to the sharp 
dependence of the quasiparticle tunneling rates on energy, and is limited by this 
sharpness. Note that unlike a normal SET or a tunnel junction NEMS, the scale of 
T e ff here is not simply set by the drain-source voltage. 



4. Back-action of incoherent Cooper-Pair tunneling 

4-1. JQP Process 

We now examine the back-action of incoherent Cooper-pair tunneling processes in a 
SSET. The simplest such process is the Josephson quasiparticle resonance (JQP), a 
transport cycle which has been studied extensively in the context of superconducting 
qubits [TBJ [T7J EU l3*2*j . The first step of a JQP cycle involves a Cooper-pair tunneling 

Note that environmental voltage fluctuations in the SET will have a similar effect of smearing out 
the energy dependence of the quasiparticle tunnel rates, and can also be parameterized by tqp. 
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Figure 1. a) Schematic of the JQP process, showing each of the three charge transfer 
steps, b) Schematic of the DJQP process. 
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through one junction of the SET, changing the charge state of the SET island by 
two electrons. This is followed by two quasiparticle tunneling events through the 
other junction which return the SET island to its initial charge state; the cycle then 
repeats (see Fig. UK) . Unlike sequential tunneling (of either electrons or quasiparticles), 
incoherent Cooper-pair tunneling is a quantum resonant process; as such, its back-action 
will have unique signatures. 

For definiteness, we focus on a JQP resonance where a Cooper pair first tunnels 
onto the island from the left, taking the island charge from n = to n = 2. We then 
have a quasi-particle tunneling event involving the right lead, taking the island charge 
from n = 2 to n = 1 (described by a rate r a ). Finally, another quasiparticle tunneling 
occurs in the right lead, taking n = 1 to n = (described by a rate T^). There are 
three important energy scales here: 5, the detuning of the Cooper-pair resonance, and 
E a (Eb), the amount of energy driving the quasi-particle event r a (r&). In the case of 
symmetric SSET junctions, these energies are given simply by: 

5 = Ef irua - E mitial = AE c (l - A/") - eV DS (25a) 
E a = 2£ c (3/2 - A/") + eVns/2 = 5/2 + eV DS + E c (25b) 
E b = 2E c {\/2 - AO + eV D s/2 = 5/2 + eV DS - E c (25c) 

The modification for non-symmetric junctions is straightforward, see e.g. [33]. The 
quasiparticle tunneling rates T a ,Tb are determined from E a , Eb and Eq. (|2"2"j). Note 
that as E a > Eb, one always has r a > r&. 

To calculate the current near the JQP feature, the standard approach is to calculate 
the equation of motion for a reduced density matrix describing the SSET island charge; 
this is done to lowest non-vanishing order in perturbation theory in tunneling in the 
junctions ^Bj- The approximation here is that the tunneling is weak, which requires the 
dimensionless conductance g of the junctions to be much smaller than 2tc (we consider 
the case where both junctions have equal conductances for simplicity). This density- 
matrix approach can be extended to calculate the quantum charge noise of the island, as 
discussed in Ref . j2H] • One again uses the "qubits as spectrometer" idea discussed in the 
previous section, and studies a system where the SSET island is coupled to a two level 
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system (TLS). In complete analogy to the calculation for the charge noise of sequential 
tunneling j2H|, the weak-coupling, long time dynamics of this system, evaluated for 
different TLS splitting frequencies, yields the quantum noise spectrum of n. 

In what follows, we present results for the transport and noise properties for the 
JQP cycle obtained from this approach. Three main interesting effects are found: i) 
JQP can lead to an extremely low effective temperature, meaning that it can be used 
for active cooling; ii) JQP can lead to negative damping effects; the resulting instability 
leads to an interesting regime of strong electro-mechanical feedback (to be described in 
detail in Sec. 5); iii) JQP can be used to reach the quantum-limit of position detection 
within a factor of two. 

The current for the JQP process is given by [Ifi] : 

/[AT, V ds ] = 2e ^fVo-lt /TM ( 26 ) 

A6 2 + {hT a y + Ej(2 + T a /T b ) 

Here, Ej = gA/8 is the Josephson energy which sets the strength of coherent Cooper- 
pair tunneling. Note that both T a and Tj, in the above equation depend only weakly 
on M and Vds near the JQP resonance; the main dependence of the current on these 
parameters is through 5, the energy detuning of the Cooper-pair resonance. To obtain 
a heuristic understanding of the JQP cycle, it is useful to consider Eq. (|26|) in the 
limit of small Ej/(KT a ); this limit is often approached in real experiments. In this limit, 
Cooper-pair tunneling is the limiting step of the JQP cycle, and one finds the simple 
result that / is 2e times the incoherent Cooper-pair tunneling rate from the charge state 
n = to n = 2: 

J^2 e xr cfTW =2ex(fi) 2 y2 + ( ^ /2); (27) 

The incoherent Cooper-pair tunneling rate has the usual form expected from Fermi's 
Golden rule: a transition matrix element squared (i.e. (Ej/2) 2 ) times a density of 
states. Here the density of states is that of the final n = 2 charge state, which is 
lifetime-broadened by the quasiparticle transition T a . 

We turn now to the back-action of the JQP process. For simplicity, we will focus on 
operating points near the center of the JQP resonance, and make use of the smallness of 
the dimensionless conductance g. In this regime, one can safely neglect the quasiparticle 
contributions to back-action damping (i.e. 9^r a , c^T;,), as they are higher order in 
g. The dominant contribution to back-action damping will instead arise from the 
oscillator's ability to modify the condition for resonance; as we will discuss, this is in 
complete analogy to damping in a number of different quantum resonant systems. The 
neglect of quasiparticle damping will allow us to derive some simple analytic expressions 
characterizing the back-action. 

In the usual limit of a small oscillator frequency (i.e. Q <C r a , T&, Ej), the back- 
action of an SSET biased near the JQP resonance is described by: 

(hT a ) 2 + 45 2 

B e// = 165 ( ' 
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A5 2 + c 1 (hT a ) 2 + c 2 E 2 



A 2 " T 2 b [(45 2 + (hT a ) 2 + c 3 E 2 ^ 



J j) 



-rnq _ T a A5 2 + c 1 (hT a ) 2 + c 2 E 2 

Tresp ~ Ak T b (T a + 2T b )' A5 2 + (hT a y + c,E 2 [ ^ C) 

where the dimensionless coefficients ci,c 2 and C3 are determined by the asymmetry 
between the two quasiparticle tunnel rates: 

ci = 1 + 4^ + 8 (^) 2 (29a) 



r a vr a/ 

c 2 = l + 4^ + 4(^) 2 (296) 

c 3 = 2 + ^ (29c) 
1 b 

Note that the above expressions are valid for arbitrary ratios of Ej/T a and T a /Tb. 

Several comments are in order. First, note that whenever the detuning of the 
Cooper-pair resonance is negative (i.e. 8 < 0), both T e ff and 5 become negative. This 
behaviour has a simple physical interpretation [2H]- From Eq. (|25a|) . we see that 5 > 
means that Vds is smaller than what is needed to have the Cooper-pairs on resonance. 
In this case, if the tunneling Cooper-pairs can absorb energy from the oscillator they 
can move closer to resonance, while if they emit energy to the oscillator, they will move 
further from resonance. The net result is that the SSET prefers to absorb energy from 
the oscillator, which implies both 7 > 0, T e ff > 0; this is the same behaviour exhibited 
by an equilibrium bath. In contrast, when 5 < 0, the situation is reversed. Vds is now 
larger than what is needed to have the Cooper pairs on resonance, and the SSET prefers 
to emit energy to the oscillator in order to move the Cooper pairs closer to resonance. 
This is the meaning of negative damping and negative temperature: our effective bath 
(i.e. the SSET) prefers to excite the oscillator rather than to absorb energy from it. 
Using Eq. (fTUj) . we see that the effect of negative damping is to reduce the total 
damping compared to that provided by the equilibrium bath, and increase the oscillator 
temperature. For 7 + 7 < 0, the total damping of the oscillator becomes negative, and 
we enter an interesting unstable regime. The properties of this regime are the subject 
of Sec. 5. 

The fact that the sign of 7 and T e ff change as we tune the detuning 5 through the 
resonance is a generic feature of the noise properties of a number of different resonant 
systems. For example, identical effects occur in an optical Fabrey-Perot cavity where 
one mirror is flexible fOlEHlE]- If such a cavity is driven by laser light which is slightly 
detuned from the frequency of the cavity, the resulting "bath" of cavity photons will both 
damp and heat the mirror. The detuning parameter in this Case is 8 — Tl(lis 'cavity ^laser)] 
when 8 < 0, photons from the laser can become resonant with the cavity if the y give 
up energy to the mirror, and one gets negative damping. This analogy can in fact be 
taken further: the effective temperature for the optical cavity system is identical to the 
expression found for JQP, c.f. Eq. ()28aj) [2*Tj . For the JQP process, T e ff is set by the 
detuning 8 and the quasiparticle rate r a ; T a is the inverse lifetime of the resonant n = 2 
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charge state. The expression for the optical cavity system is identical, with T a being 
replaced by the inverse lifetime of a photon in the cavity [21]. The correspondence is a 
direct consequence of the fact that in both cases, the motion of a mechanical degree of 
freedom modifies a resonance condition. 

While the analogy to the optical cavity system is certainly useful, there are some 
respects in which it is different from JQP case. Consider the expression for the response 
time r resp of our effective bath. As expected, r resp scales as the lifetime of the resonant 
state (i.e. 1/T). More surprising is the fact that r resp depends on the detuning 5: our 
effective bath responds slower when we are on resonance. This is in sharp contrast to 
what is found for the cavity plus mirror system discussed above; there, the response 
time is always set by the cavity ring-down time, regardless of detuning [T§1 HU[ l2*Tj. 

Finally, we wish to emphasize the smallness of T e ff for the JQP process: unlike 
the normal metal SET or tunnel junction, the effective temperature here is not set by 
the source-drain voltage, but rather by a much smaller scale, T a . As with the case of 
sequential quasiparticle tunneling, we see that the T e ff of a mesoscopic conductor is by 
no means always set by the drain source voltage. The effective temperature is rather 
a measure of the asymmetry between energy absorption and emission. In the present 
case, the smallness of T a means that incoherent Cooper-pair tunneling can be used for 
substantial cooling of the resonator when 5 > 0; we give numerical estimates of the 
magnitude of this cooling effect in Sec. 4.3. 

4-2. DJQP back-action 

We now turn to the back-action of a slightly more complex incoherent Cooper-pair 
tunneling process, the Double Josephson Quasiparticle resonance (DJQP). This cycle 
consists of two incoherent Cooper-pair tunneling events in series, one in each junction 
(see Fig. [TJd) ; its back-action properties (in the context of qubit measurements) were 
studied in [2E]- The DJQP process is of interest as it occurs at a smaller drain-source 
voltage in the SET than the JQP process; as we will see, it also has more pronounced 
back-action, and can have a greater cooling effect. 

For definiteness, we consider the following cycle: a) Cooper-pair tunnels from left, 
taking the island charge from n = to n = 2; b)Quasi-particle tunnels from island to 
right lead, taking n — 2 to n = 1; c)A Cooper-pair tunnels from the island to the right 
lead, taking n = 1 to n = -1 ;d)A quasiparticle tunnels onto the island from the left 
lead, taking n = — 1 to n = . At this stage, the cycle repeats. Note that there are two 
Cooper-pair tunneling events (one in each junction), and two quasi-particle tunneling 
events (one in each junction). We denote the first quasiparticle transition rate (n = 2 
to n — 1) by r a , and the second by r^,. We also denote the first Cooper-pair transition 
(n = to n = 2) as "A", and the second transition {n = 1 to n = — 1) as "B". As with 
the JQP process, we have to worry about the following volt age- dependent energies: 



5 A =AE c {l-M)-eV DS 
5 B = - 4E C (0 - Af) - eV DS 



(30a) 
(306) 
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E a = 2E c {?>/2 -Af) + eV DS /2 (30c) 
E b = - 2E c (-l/2 -Af) + eV DS /2 (30d) 

Here, 5a,o~b are the detunings of the two Cooper-pair transitions, and E a , E b are the 
energies driving the two quasiparticle transitions. The rates T a , T b of these quasiparticle 
transitions are related to the corresponding driving energies by Eq. f!22|) as before. Note 
that the center of the DJQP resonance (where both 6 a arid 8b are 0) occurs when both 
M = 1/2 and eV ds = 2E C . 

An density-matrix approach identical to that used for the JQP may be used to 
describe the transport and noise properties of the DJQP. While this technique is valid 
for arbitrary ratios Pb/r a and Ej/T a , we will make some additional mild approximations 
to obtain results that may be easily interpreted. First, we again consider a SSET with 
identical junctions, and note that at the center of the DJQP resonance, T a = T b . As 
we are interested in the behaviour near the resonance center, we will take r a = = T 
throughout. In addition, we will consider the limit Ej -C T; this limit is approximately 
realized in many experiments, and leads to a great simplification in the form of the 
resulting equations. 

In this limit, the Cooper-pair tunneling is the rate limiting step in the DJQP cycle; 
we find that the average current is given by 3e times the series addition of the incoherent 
Cooper-pair tunneling rates for the "A" and "B" transitions: 



3e x 



1 



-i 



+ 



3e 



EjT 



2 (nvy 



2S\ + 26% 



(31) 



\Tcpt(8a) Fcpt(o~b) 
The incoherent Cooper-pair tunneling rate Tcpt{5) is defined in Eq. (12 7j) . 

In calculating the back-action for the DJQP process, we make again make 
the assumption that we can ignore the quasiparticle contribution to damping; as 
discussed, near resonance, this contribution leads to terms which are higher order in the 
dimensionless conductance g. For the effective temperature, we again find a remarkably 
simple result: 

i r 1 

+ 



kE.Teff 



k R T, 



BJ-eff,B 



+ 



165 



B 



ksT e ff,A 
IGS A 

(hTf + A8\ ' (hT) 2 + 45 2 B 
[(hTf + A5\] [{hY f + 4§] 

lQ[6 A + 6B}[(hT) 2 + 45 A 6B} 1 ' 

We see that the effective temperature for the DJQP process (in the small Ej limit) 
is simply given by the series addition of the effective temperature for each individual 
Cooper-pair resonance; the scale for the minimum effective temperature is again set by 
T (i.e. the lifetime of each of the two resonant states). The equation for the back-action 
damping 7 is slightly more complicated: 



777,7 
^2" 



) Sa+_5b_ 



((ftr) 2 + A5\){{hT) 2 + A5 2 B ){{hT) 2 + ASaSb) 
[{hT) 2 + 25 2 A + 25 2 B f 



(33) 
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Similar to the situation with JQP, the effective temperature (and damping) can 
become negative if the net preference of the two resonances is to emit energy to the 
oscillator (i.e. the series addition of T e ff t A and Te//,B must be negative). To see exactly 
where negative damping and temperature occur in terms of SSET operating points (i.e. 
Vds and A/"), it is useful to re- write the expression for T e ff. We define: 

AE C M = 4E c {l/2 + 5M) = 2E C + A M (34a) 
eV DS = 2E C + eSV DS = 2E C + Ay (345) 

Ajv and Ay represent the energy detuning that result when we move, respectively, Af 
or Vds from the center of the resonance. With this notation, we have: 

t T 1 [r 2 + 4(Ay + A^) 2 ] [r 2 + 4(Ay - A^) 2 ] 

kBleff - -32 Ay ' T 2 + 4 [(Ay) 2 - (A^) 2 ] {6b} 

For small detunings from the center of the resonance (Ay, A^ -C T), we see that the 
sign off the effective temperature is completely determined by the value of the drain- 
source voltage. If Vds is set to be below the center of the DJQP resonance, Ay is 
negative, and hence both T e ff and the damping 7 are positive. In this case, the net 
preference of the two resonances is to absorb energy from the oscillator. In contrast, 
if the drain-source voltage is tuned to be higher than resonance, we get both negative 
temperature and negative damping. Now, the DJQP process can be brought closer to 
resonance by emitting energy to the oscillator. Finally, note if the drain source voltage is 
exactly at the center of the resonance (i.e. Ay = 0), T e ff tends to infinity, meaning that 
there is no damping, and no asymmetry between absorption and emission; this holds 
regardless of how much one moves the gate voltage M off-resonance. In this special case, 
one Cooper-pair transition always prefers to absorb energy, the other to emit; the net 
result is no preference between absorption and emission. In Fig. 121 we show a contour 
plot of T e fj for DJQP using typical device parameters, and compare it against a plot of 
the current. 

It is instructive to compare the magnitude of detector-induced damping for the JQP 
and DJQP processes in the small Ej limit we consider. In both cases, the minimum value 
of \T e ff\ is oc T, where T is the quasiparticle tunneling rate. However, the corresponding 
damping will be much larger in magnitude for DJQP versus JQP: comparing Eqs. (|286|) 
and (}33|) . we see that the back-action damping in the DJQP case is enhanced over that 
at JQP by a large factor (T/Ejf. The reason for this is simple to understand, once we 
recall that 7 oc S n (0)/T e ff, where S n (0) is the zero- frequency charge noise of the SSET 
island. For both JQP and DJQP, this noise has a telegraph-noise form in the small Ej 
limit. For JQP, we have effective telegraph noise between the n = and n = 2 charge 
states, where one rate (the quasiparticle transition rate T) is much larger than the other 
(the incoherent Cooper-pair tunneling rate T C pt)- As a result, S n (0) oc T C pt/F 2 . In 
contrast, for DJQP, we have effective telegraph noise between the n = and n = 1 
charge states; the two telegraph rates are both incoherent Cooper-pair tunneling rates 
(i.e. the "A" and "B" transitions), and are roughly equal. In this case S n (0) oc 1/T C pt- 
Thus, as both processes have similar effective temperatures, the enhanced charge noise 
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Figure 2. a) Contour plot (with equally spaced contours) of current near DJQP 
resonance, as a function of gate voltage detuning (A/V = TV — 1/2) and drain-source 
voltage detuning (AVds — Vds ~ 2Ec/e) from the resonance center. We have chosen 
typical SSET parameters A = E c = 0.2meV and g = 0.2. b) Contour plot of (T osc )~ , 
the inverse of the oscillator's temperature (c.f. Eq. iJHJ}), for the same range of AVds 
and AAf. We have taken typical oscillator parameters Q,/(2-k) = 27MHz, Qq = 10 4 
(i.e. intrinsic oscillator quality factor), d — 300nm (i.e. oscillator-SET island distance), 
C osc — O.lCs and V osc = 10V (c.f. Eq. (fTTjl l: we have also chosen a bath temperature 
of T = 500mK. Yellow indicates unstable regions where the oscillator temperature is 
negative. In the regions of cooling (blue), the lowest achievable oscillator temperature 
is less than 50mK, a factor of ten lower than the oscillator's temperature at zero 
coupling to the SET. 



of DJQP explains its enhanced damping. This enhancement of back-action damping 
should make it easier to see back-action effects of DJQP in an experiment. 

Finally, consider the response time r resp of the effective bath presented by DJQP: 
it is given by: 

= ^cpt{&a) + Fcpt{$b) 

7~resp 

2Ej {hTf [{hTf + gj| + 26%] 
" h 2 Y [{hVY + A5\]{{hTf + A5l} 1 ] 

The average response rate, l/r rep , is simply given by the sum of the rates of the two 
Cooper-pair tunneling events in the cycle. In the small Ej limit we are considering, 
this timescale is much longer than what was found for JQP (c.f. Eq. (J28cj) ). where one 
always had r resp ~ ^. 

As a final note, we should remark that none of the approximations used here (small 
Ej, no quasiparticle contribution to damping) are necessary to the method used. We 
have made these approximations only so that we could give a clear heuristic picture of 
the physics underlying the back-action effects found. 
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4-3. Typical Cooling Temperatures 

Eq. (jlOj) tells us that if the back-action damping dominates the intrinsic, non-back- 
action damping of the oscillator (i.e. 7 ^> 70), then it should be possible to cool the 
oscillator to the effective temperature of the SET, even if this is much lower than T 
(i.e. the temperature of the oscillator when it is decoupled from the SET). To achieve 
7 ^> 7o, one should start with an oscillator with a high intrinsic quality factor, and 
engineer structures where the SET-oscillator coupling is relatively strong (e.g.. have the 
SET island be in close proximity to the oscillator) Note that this condition does not 
imply any violation of the weak-coupling assumptions of our treatment: it is entirely 
possible to have the coupling A be weak-enough that linear response is still valid, yet 
strong enough that 7 ^> 70 . 

To give an idea of the kind of cooling that should be achievable using incoherent 
Cooper-pair tunneling, it useful to give some numerical estimates using typical SET 
parameters. We assume Eq = A = 0.2meV, in agreement with the devices used in Ref. 
[31]. For the JQP, we saw that the minimum possible positive value of T e ff was given 
by T a /4. Choosing Vns = 2A + Eq (the threshold voltage for the JQP process), we find 
for these parameters: 

g x 350mK (37) 

Thus, for g = 0.2 (i.e. a junction resistance of !30kQ ), the minimum effective 
temperature is approximately 70mK; again, this is the minimum temperature we could 
cool our oscillator to using the JQP process. Note that this T e ff is much lower than the 
scale set by the voltage: cVbsI^b — 7. OK. 

Turning to the DJQP, and now assuming the typical condition Ej < T, we see from 
Eq. (J3*2"|) that the minimum possible positive value of T e ff is T/8, a full factor of two 
smaller than for the JQP process. If we again take Ec = A = 0.2meV, we find: 

g x 125mK (38) 

For systems with Ej < T, cooling should be much easier to achieve using DJQP. Not 
only is the T e ff lower, but, as already discussed, the back-action damping due to DJQP 
is larger than that of JQP by the large factor (T/Ej) A . 

As a final caveat, the above estimates assume that the only source of broadening of 
the Cooper-pair tunneling resonances is quasiparticle tunneling. In reality, additional 
effects could further broaden these resonances and increase the minimum achievable 
temperature. Perhaps the most important of these will be voltage fluctuations associated 
with the environmental impedance seen by the SET island. These effects will ultimately 
limit how much one can lower T e ff by simply lowering g. 



T eff 
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4-4- Quantum limited measurement with incoherent CPT 

There is an additional important point that needs to be made about incoherent Cooper- 
pair tunneling processes: unlike sequential tunneling in a normal SET, these processes 
may be used to do near quantum-limited displacement detection [2BJ- As discussed 
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Figure 3. Noise temperature for the JQP and DJQP processes (scaled by Ml/ 2) versus 
Ej/T, where Ej is the Josephson energy, and T is the quasiparticle tunneling rate. 
We have assumed an optimal choice of coupling strength A, identical SET junctions, 
and that both quasiparticle tunneling rates in each cycle are the same. 



extensively in j22] , reaching the quantum limit requires a detector which has a minimal 
amount of back-action noise relative to the amount of information it provides. The 
measure of this relative back-action is the noise temperature Tjy of the detector, which 
gives a measure of how much noise is added to the signal by the detector. The quantum 
limit is that /cbT/v cannot be any smaller than hQ/2, where Q is the oscillator's frequency 
We are interested in the usual limit where the detector noise correlators are 
frequency independent on frequency scales relevant to the oscillator (i.e. Q, 7). After 
optimizing the strength of the coupling A (to balance back-action and intrinsic shot 
noise contributions to the total noise), one has 



k B T N 

X = 



> , (39 ) 



IX) 



htt/2 ^ (hdi/d. 

Here, Sp = A 2 S n is the back-action force noise spectrum of the detector, Sj is the 
zero frequency current noise of the SSET, and dl/dx = -{dl/dM)^ is the gain of our 
detector: how strongly does the average current respond to changes in the oscillator's 
position. 

Thus, achieving the quantum-limit of displacement detection requires that x 
approach one + . For a normal-metal SET in the sequential tunneling regime, x is 

+ Note that if our detector was used to do a weak QND measurement of a qubit, 1/x 2 represents 
the ratio between the rate at which information is acquired by the measured versus the back-action 
dephasing rate. 
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proportional to 1/g 2 , where g is the dimensionless conductance of the junctions. By 
assumption, g is a small parameter- this is what allows us to use perturbation theory in 
the tunneling! One thus has that sequential tunneling in a normal SET is far from being 
at the quantum limit; there is considerable excess back-action noise, far beyond what 
is required by quantum mechanics. In contrast, for JQP and DJQP, \ is order 1, and 
is independent of g [2E|. Using the same density matrix approach used to calculate the 
charge noise, one can calculate the zero frequency current noise of both JQP and DJQP 
[S21 12H1 EE! • Shown in Fig. El is the reduced noise temperature x f° r JQP an d DJQP, 
evaluated as a function of the ratio Ej/Y, and for a bias point which yields a maximum 
gain. Note that the ratio Ej /Y can in principle be controlled if one can control the ratio 
of Eq to A when fabricating the SET. We see when that Ej ~ Y, both DJQP and JQP 
come close to reaching the quantum limit, with the DJQP process doing slightly better. 
In the more usual case where Ej <^Y, the JQP process is much better at approaching 
the quantum limit than the DJQP process; this is because of the relative suppression of 
charge noise associated with JQP when Ej C T. 

5. Negative Damping and Strong Electro-mechanical Feedback 

In this last section, we describe in more detail the interesting negative-damping 
instability that can be brought about by incoherent Cooper-pair tunneling in a SET 
plus mechanical resonator system. As discussed, the negative back-action damping 
associated with either the JQP or DJQP transport cycles can, for a sufficiently strong 
coupling, drive both the total oscillator damping and its temperature T osc (as given 
by Eq. ()10j)) below zero (e.g. see the yellow regions in Fig. Efc>). In this regime, 
the SET will continually dump energy into the oscillator, causing the amplitude of 
the oscillator's motion to continually increase. Eventually, this amplitude will become 
so large that there will be an effective strong coupling between the oscillator and the 
SET- the motion of the oscillator will strongly effect the dynamics of the SET. In this 
regime, the weak coupling approximations we have been making clearly break down. 
We nonetheless wish to describe the properties of the eventual stationary state achieved 
in this regime. Note such negative-damping instabilities are also well known in the 
analogous optical cavity system [371 EHJ EH| ■ 

To deal with the feedback occurring in this effective strong coupling regime, we 
will make the usual simplifying assumption that the oscillator is much slower than the 
SET: O C T, Ycpt'i this condition is usually more than satisfied in experiment. We can 
then make use of this separation of timescales in manner somewhat analogous to the 
usual Born-Oppenheimer approximation. As in our original weak-coupling description, 
at each instant in time the SET acts as an effective thermal bath on the oscillator, 
and is characterized by a fluctuating back-action force with spectral density S$f and a 
damping rate 7; the oscillator is thus still described by the Langevin equation Eq. 0. 
Now, however, the parameters of the bath will depend on the instantaneous position 
of the oscillator: this is how we deal with the strong-coupling aspect of the problem. 
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Thus, the spectral density of the back-action force 5 fit) appearing in Eq. (jHJ) will be 
x-dependent (i.e. S$f = Sgf[x)), as will the damping coefficient 7 (i.e. 7 = j[x]). The x 
dependence of both these quantities arises completely through their dependence on the 
SET dimensionless gate voltage A/"; we can absorb the coupling Hamlitonian of Eq. (0) 
into the SET charging energy Hamiltonian (c.f. Eq. (|15jl) by writing: 

N[x] = m + 7^ x = — + 4r x - ( 4 °) 

Again, the key assumption here is that the SET response to changes in x is much faster 
than the evolution of x itself. 

As we show in what follows, the net result of this approach is that we may describe 
our coupled system with a classical Fokker-Planck equation where both the diffusion 
and damping constants have an explicit x dependence. Heuristically, for large oscillator 
amplitude, the Cooper-pair detuning S and hence the back-action damping will depend 
on the position of the oscillator; in the unstable regime of interest, the stationary state 
of the oscillator will correspond to an amplitude of oscillation so large that the oscillator 
experiences equal amounts of positive and negative damping during one period of its 
motion. More quantitatively, we find that the stationary state of oscillator is naturally 
characterized by an energy- dependent temperature; for a variety of regimes, the energy 
distribution of the oscillator is a gaussian, and hence highly non-thermal. Note that a 
similar approach to potential strong-feedback behaviour was used in Ref . ^T] ; however, 
the particular NEMS systems studied in that work did not exhibit any negative damping, 
and thus there was no effective strong-coupling regime. 

Our starting point is thus the Langevin equation of Eq. Q with x dependent 
damping and back-action force terms. As we are interested in the limit where the 
oscillator is much slower than the source of the back-action noise, we can treat the 
latter as being white. We may then convert our Langevin equation to a Fokker-Planck 
equation for the oscillator's phase-space density w(x,p;t) in the usual way [40; we 
obtain: 



d 

di W 



p d d 



m dx dp 



+ — (mtt 2 x + (70 + j[x])pf) 



q2 

w + (D Q + D[x))—w. 



(41) 

Here 70 is the intrinsic (SET-independent) damping of the oscillator, D = mjokBTo 
describes momentum diffusion due to the equilibrium bath, and the x-dependent 
diffusion constant D[x] is determined from back-action force spectrum: 

D[x] = -S sf [uj = 0; x] = m ■ 7 [a;] • k B T eff [x] (42) 



In the case of JQP, 7(2) and T e ff[x] are given respectively by Eqs. (j28fc|) and 
with the substitution Af — ► M[x\\ for DJQP, one would use Eqs. (jHBj) and (J32j) . Note 
that there is no ambiguity in interpreting Eq. (JHJ) with an x-dependent back-action noise 
spectrum: both the Ito and Stratonovitch interpretations of this stochastic differential 
equation yield the same Fokker-Planck equation. Also note that we are neglecting the 
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x-dependence of the conservative part of the average back-action force; this is justified 
in the experimentally relevant limit of a high intrinsic quality factor Q . In this limit, 
one can have the magnitude of back-action damping be comparable or greater than 
the intrinsic damping of the oscillator, while at the same time have the back-action 
spring-constant modification Ak be much much smaller than mQ 2 . 

To make further progress, we consider the experimentally-relevant weak-damping 
limit of Eq. (JIT]) (i.e. 7 , <C fi), and follow the approach of Kramers [JT] (and 

more recently, Blanter et. al [TT]). We first re-write x and p in terms of the oscillator's 



energy E and phase 9 (x = y2E/ksin9, p = \ / 2mE cos 9), and then convert Eq. (|41|) 
into an equation for w(E,9). Working to lowest order in the dissipative terms (both 
damping and diffusion), this equation can then be recast as an equation for the oscillator 
energy distribution w(E) = J 2n d9 w(E, 9); one uses the fact that in the absence of any 
dissipation wq(E,9) = wo(E)/2n. After some algebra, this procedure yields: 

d x d „ / /rnN D + D(E) d \ , 

-w(E; t) = -E U + 1( E) + «,<*; t) (43) 



The energy-dependent back-action damping and diffusion constants here are defined as: 

cos 2 9 (44) 



f 2w d9 
D[E] — 2 [ 2 * ^-D 



x = J2E/ksm9 



x = J2E/ksm9 



■ cos 2 9 (45) 



/o 2tt 

As could have been expected, they are given by averaging ~f[x] and D[x] over the phase 
of the oscillator at fixed amplitude; what is perhaps more surprising is that one must 
take a weighted average, with a weighting factor which is proportional to p 2 . Note 
that at E = 0, j(E) and D[E] coincide respectively with the back-action damping and 
diffusion constant in the weak-coupling (i.e. zero feedback) theory. 

It is now trivial to solve for the stationary energy distribution of the oscillator. 
From Eq. (JIB]), we have: 

/ r E dE' \ 
w(E) stat — N exp I — / — (46) 



k B T osc [E'} 

where is a normalization constant, and the energy-dependent effective temperature 
T osc (E) is given quite naturally by: 

ksf^E) = D , 0+ fW (47) 
m(7o + 7[S]) 

We see that the stationary state is characterized by a generalized Boltzmann distribution 
having an energy- dependent effective temperature. In the weak-coupling, no negative 
damping case, we can neglect the energy dependence of kBT[E], and Eq. coincides 
with a thermal distribution and with the results of the weak-coupling theory (c.f. Eq; 
(JTUJ)). In the strong-feedback regime of interest, the total damping ( i y + r y(E)) is negative 
at E = 0; as we increase E, the motion of the oscillator smears out the back-action 
contribution to the damping, and for very large E, the total damping tends to 7 > 0. 
There will thus be a critical E = Eq where (70 + 7(^0) ) w iU P ass through zero. It follows 
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Figure 4. Energy-dependent back-action damping of the JQP process, c.f. Eq. I|44|l . 
plotted as a function of energy, and using SSET parameters identical to those used in 
Fig. 2. The characteristic scale for feedback effects, E* , is given in Eq. I|5U|) . Note that 
the chosen scaling makes this plot independent of the coupling strength A. The red 
curve and blue curves correspond to two different Cooper-pair detunings 8 (i.e. SSET 
operating points) as labeled; note that j(E) can have a non-monotonic dependence on 
energy. The inset shows the dependence of the back-action damping at zero oscillator 
energy (c.f. Eq. I|286|l ) on the Cooper-pair detuning. 



directly from Eq. (}46|) that w(E) will have a maximum at Eq; moreover, in the vicinity 
of E , w(E) will look gaussian: 

w(E) stat ~ iVexp (- {E ~_y 2 ) (4f 



a 



Dq + D[E 
2m 



2a 2 

d'y 



dE 



E=E . 



-1 



(49) 



We thus find that in the strong-feedback regime, the oscillator's energy distribution 
function has a maximum at a non-zero energy, in sharp contrast to an equilibrium 
distribution. The energy at which the maximum occurs corresponds, as expected, to an 
oscillator amplitude large enough that the oscillator experiences zero average damping 
during each period of its motion. The width of the distribution near this maximum is set 
by both the back-action noise and the energy-sensitivity of the damping via Eq. ()49|) . 

We now apply the above results, which are quite general, to the specific instability 
brought on by incoherent Cooper-pair tunneling. We first consider the case of 
an oscillator-SSET system operated in the negative damping regime near the JQP 
resonance; as usual we consider the limit of small dimensionless junction conductances, 
and neglect the quasiparticle contribution to the back-action damping. In this limit, all 
the feedback effects of the oscillator on the SET will be due to the dependence of the 




Figure 5. Effects of JQP back-action in the negative-damping regime. We have taken 
a bath temperature To = 50 mK, an intrinsic oscillator quality factor Qo = 10 5 , and 
a SSET operating point near the JQP resonance which maximizes the effect of the 
back-action damping (S = — 0.25?z.r a ). All other parameters are the same as in Fig. 2. 
Left panel: Average energy of the oscillator E avg as a function of the coupling voltage 
V osc ; the dashed curve is the predictions of the weak-coupling theory, c.f. Eq. I jlOj l. 
while the red curve is from the strong-feedback theory. The total oscillator damping 
becomes negative at V osc = 7.1 V. Right panel: Oscillator energy distribution w(E) as 
determined by Eq. . As V osc becomes large enough for the strong-feedback effect 
to take hold, the energy distribution becomes gaussian-like. 



Cooper-pair detuning 5 (c.f. Eq. (|25ajl ) on the oscillator position x. The characteristic 
energy scale E* for feedback effects will correspond an oscillator amplitude large enough 
that the corresponding oscillations in 5 are equal to the resonance width, hT a . From 
Eqs. (HH and (pZ5a|) , one has: 




E* sets the characteristic scale for variations in ^[E] and -D[-E]- Shown in the inset of 
Fig. 4 is the back-action damping 7 at zero energy (i.e. ignoring feedback effects), c.f. 
Eq. (|286Jl . as a function of 5. The main plot shows the energy dependence (scaled by 
E*) of the the damping for two choices of Sq = 5[x = 0]. Note that j(E) can have a 
non-monotonic dependence on E. 

While measuring the full distribution w(E) of the oscillator's energy may be quite 
difficult, the average oscillator energy E avg can be obtained directly from experiment. In 
the experiment of LaHaye et. al [2j, this quantity was obtained by extracting (x 2 ) from 
the current noise of the SET, and then invoking the equipartition theorem; a similar 
approach could be used in the strong-feedback regime discussed here *. In Fig. [EI a), we 
plot E avg as a function of the coupling voltage V osc ; we use typical device parameters 
(see figure caption), and choose 5 = —0.25hT a to maximize \j(E = 0)| (i.e. same as the 
blue curve in Fig. 4.). For small V osc , the intrinsic damping of the oscillator dominates 

* As the stationary state is characterized by a flat distribution of the oscillator phase 8, one again has 
the relation E avg /k — (x 2 ), as would hold in equilibrium 




Figure 6. Effects of JQP back-action in the negative-damping regime. We use the 
same parameters as FigEJ but now choose an SSET operating point near the JQP 
resonance which yields a non-monotonic j(E). (S = —TiT a ). Left panel: Average 
energy of the oscillator E avg as a function of the coupling voltage V osc , labeled as in 
FigO Note that the strong-feedback effect sets at V osc much lower than what is needed 
to make the weak coupling theory diverge. Right panel: Oscillator energy distribution 
w{E) as determined by Eq. l(I3Jl. 



the back-action damping, and there is no strong feedback effect. The oscillator is in 
a thermal state, with a temperature T osc given by Eq. ([TO]), and E avg = k B T osc . As 
we increase the coupling voltage, the back-action damping increases in magnitude (i.e. 
becomes more negative), and E avg = T osc correspondingly increases. At a critical voltage 
V osc = 7.1 volts, the total damping becomes negative and the weak-coupling expression 
for T osc diverges. The strong-feedback effect takes hold here, and we see that average 
energy of the oscillator remains finite. There is also a corresponding change in the 
shape of the energy distributions, as can be seen in Fig. E3b). Note that for a large 
range of V osc , E avg remains approximately constant. This is the result of two competing 
tendencies: increasing V osc increases the overall magnitude of the back-action damping, 
but also increases the sensitivity of the Cooper-pair detuning 5 to E. 

In Fig. we plot the same quantities as0 but now for 5 = —TiT a . In this case, 
the non-monotonic nature of j(E) leads to energy distributions w(E) which have two 
local maxima. Also note that in this case, the strong feedback effect sets in well before 
the weak-coupling expression for T osc diverges. 

Finally, in Fig we plot results for an oscillator SSET system operated in the 
negative-back action damping regime near the DJQP resonance. We have again chosen 
realistic parameter values (see figure caption). We see that the feedback effect is 
stronger here; as discussed, this is a consequence of the enhanced back-action damping 
of DJQP versus JQP. Note that for both the JQP and DJQP strong feedback results 
presented here, the maximum amplitude of the oscillator is still much smaller than the 
oscillator-SET gap d, and the corresponding oscillation of the Cooper-pair detuning is 
at most order T. Thus, these results do not violate our assumption of a linear coupling 
Hamiltonian, and do not involve SSET physics far from the JQP or DJQP resonance. 
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Figure 7. Effects of DJQP back-action in the negative-damping regime. We have 
taken a bath temperature To = 50mK, an intrinsic oscillator quality factor Qq = 10 4 , 
and a SSET operating point near the DJQP resonance ( AjV = 0, AV = KT/2). All 
other parameters are the same as in Fig. 2. Left panel: Average energy of the oscillator 
E avg as a function of the coupling voltage V 08C ; the dashed curve is the predictions of 
the weak-coupling theory, c.f. Eq. (|10() . The total damping becomes negative when the 
coupling voltage V osc — 1.75 volts. Right panel: Oscillator energy distribution w(E) 
as determined by Eq. (|43|l . At V osc = 1.5 volts, we still have a thermal distribution; 
increasing the coupling voltage takes one into the unstable regime and results in an 
almost gaussian distribution of energy. 



There are a number of other interesting issues related to the strong feedback regime 
discussed here. In particular, what are the dynamics of the oscillator's energy in this 
regime? As the total damping vanishes at the most probable energy Eq, one expects 
very slow dynamics. How does this then manifest itself in the output noise of the SET? 
It would also be interesting to investigate systems where the SET is not infinitely fast 
compared to the oscillator; in this case, one could expect similar kinds of multistability 
effects found in optical cavity systems jSH] • We hope to address these issues in the near 
future. 

6. Conclusions 

In this paper, we have discussed various aspects of back-action physics in SET-oscillator 
NEMS using a quantum-noise approach. In the sequential tunneling regime, we derived 
general expressions for the effective temperature, damping rate, and bath response 
time that are valid regardless of the particular system (electrons in a normal SET, 
quasiparticles in a superconducting SET, electrons in a quantum dot). In the regime 
of incoherent Cooper-pair tunneling, we described the back-action of both the JQP and 
DJQP processes, and demonstrated that both could be used for substantial oscillator 
cooling, as well as near quantum-limited measurement. We also discussed the strong 
analogy between these processes and an optical cavity with a moveable mirror. Finally, 
we discussed the regime of strong electro-mechanical feedback that can arise when the 
back-action damping becomes negative, and showed that the stationary state in this 
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regime is naturally characterized by an energy-dependent temperature. 
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